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The formalism of Siegert states to describe the resonant scattering in quantum theory is extended 
to the resonant scattering of electromagnetic waves on periodic dielectric arrays. The excitation of 
electromagnetic Siegert states by an incident wave packet and their decay is studied. The formalism 
is applied to develop a theory of coupled electromagnetic resonances arising in the electromagnetic 
scattering problem for two such arrays separated by a distance 2h (or, generally, when the physical 
properties of the scattering array depend on a real coupling parameter h). Analytic properties of 
Siegert states as functions of the coupling parameter h are established by the Regular Perturbation 
Theorem which is an extension the Kato-Rellich theorem to the present case. By means of this 
theorem, it is proved that if the scattering structure admits a bound state in the radiation continuum 
at a certain value of the coupling parameter h, then there always exist regions within the structure in 
which the near field can be amplified as much as desired by adjusting the value of h. This establishes 
a rather general mechanism to control and amplify optical nonlinear effects in periodically structured 
planar structures possessing a nonlinear dielectric susceptibility. 
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I. INTRODUCTION 



To outline the scope of the problem studied, it is helpful first to glean over the basic concepts of quantum resonant 
scattering theory. Consider the scattering problem H^ = £ "J for a quantum mechanical system described by the 
Hamiltonian H = — + V(r) where A is the Laplace operator in space, £ is the spectral parameter (the energy), 
and the scattering potential V(r) is assumed to be of a finite range, i.e., V(r) = if r = |r| > r for some r . To 
avoid excessive and unnecessary technicalities, the potential is assumed to be spherically symmetric, V(r) — V{r). 
Suppose further that only a spherically symmetric scattered wave is of interest (the s— wave). Then the scattering 
theory [2, [3] requires that in the asymptotic region r > ro an eigenfunction of the Hamiltonian H corresponding to 
an eigenvalue £ is the superposition of an incident wave "J^ ~ e lk r and an outgoing wave s : 



e ikr 



tf = tfi + tf., V S = S , k=\k\=V2£ 

r 

where S is called the scattering amplitude . The system is said to have a scattering resonance if S has a pole in the 
complex energy plane: 



£ £n ~f~ 2-1-V. 



where £ n is the resonance position (or the resonant energy) and T n is the resonance width. In this case the scattering 
cross section as a function of £ exhibits a Lorentzian profile 

a(£) ~ \Sf 



(£-£ n r+Tl 



In 1939, A. Siegert showed [3] that positions and widths of resonances for a given quantum system can be determined 
by solving the eigenvalue problem under the outgoing wave boundary condition which in the simplest case of a 
spherically symmetric scattered wave (the s— wave) reads 



H*„=£*„, (9 r (r*„)-zfcr*„) 



= 

r>ro 



Under this boundary condition the Hamiltonian is no longer hermitian, hence, the eigenvalues are generally complex, 
£ = £ n — iT n . The eigenstates \l/ n associated with eigenvalues £„ — iT n are now called Siegert states. The scattered 
wave is then proved to have the form 

n 

where the complex amplitudes a n ~ a n [^i) are homogeneous linear functionals of the incident wave and $ a is the so 
called background or potential scattering which is analytic in £ (see e.g.,0, Q). If ^ is a wave packet with a narrow 
energy distribution centered at £ — £ c (e.g., a narrow Gaussian wave packet), then the amplitude a n is significant, 
provided £ c ~ £ n . In other words, a Siegert state contributes to the scattering wave if the incident wave has a resonant 
energy £ n determined by the real part of the corresponding Siegert eigenvalue. 

Note that the scattered modes always have positive energy £ = k 2 /2 > 0. The range of £ corresponding to scattered 
modes is called the radiation continuum. Suppose that V = V(r), for simplicity. It is easy to see that the Siegert 
states contain bound states of the system £ n < and r„ = which decay exponentially in the asymptotic region 
(as k = i\J — 2£ n ) and, hence, are square integrable eigenstates of the Hamiltonian. The bound states have a discrete 
spectrum which lies below the radiation continuum (£ n < 0). If the asymptotic behavior of the potential V is relaxed 
by demanding that V(r) — > as r — ¥ oo, then there may exist bound states (i.e., square integrable solutions of the 
stationary Schrodinger equation) with positive energies £ n > 0, T n = 0. They are called bound states in the radiation 
continuum or resonances with the vanishing width (to emphasize the fact £ n > 0). Their existence was first predicted 
by von Neumann and Wigner in 1929 6]. It can further be shown that the amplitudes a n vanish if r n = 0. Thus, no 
bound state (either above or below the radiation continuum) of the system can be excited by an incident wave [l|. 

The physical significance of Siegert states with T„ > can be understood through the initial value problem for the 
time dependent Schrodinger equation id t ^ = H^ in which a wave packet initially positioned in the asymptotic 
region r > r (i.e., at the initial time t = the support of lies in the asymptotic region), propagates into the 
scattering region r < ro and passes through it giving rise to scattered waves. In the time-dependent picture, the 
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amplitude of each Siegert state (with a sufficiently small T n ) that has been excited by the incident wave packet is 
shown to decay exponentially [l|: 

*.(t) = J2 fi nW*» + *,(*) , n»(*) ~ e-^'e-r-' 

If the incident wave packet passes through the scattering region faster than the decay time t„ = l/T n of a Siegert 
state, then the outgoing wave can be observed, and it resembles to a stationary state with the energy £ n . The more 
narrow the scattering resonance is, the longer lives the corresponding Siegert state. So, Siegert states with r n < 1 
may be interpreted as quasi-stationary states of the system that can be excited by an incident wave and live long 
after the scattering process is over. 

Resonant scattering phenomena are also quite common in electromagnetic theory. In the past decade much attention 
has been devoted to experimental and theoretical studies of reflection and transmission properties of periodically 
perforated films and similar periodic planar structures in infrared and visible light. The transmission and reflection 
coefficients of such structures exhibit typical resonant peaks at a wavelength about the structure period (for a 
review see 0,0). Numerical studies of the scattering problem for an electromagnetic wave packet impinging a periodic 
structure show that a portion of the electromagnetic energyof the wave packet is trapped by the structure and remains 
in it long after the wave packet has passed the structure [9l-[l3|. Furthermore, each "trapped" electromagnetic mode 
decays slowly by emitting a monochromatic radiation whose wave length is close to the resonant wavelength. Broad 
and narrow resonances correspond to short-lived and long-lived "trapped" modes, respectively. It has been shown 
that "trapped" modes may occur either due to the structure geometry or the dispersive properties of the structure 
material 0, [ll[ . 

If two planar periodic structures are separated by a distance 2h, then the resonance positions and their widths 
depend on h so that under some conditions, there are critical values h = hb at which the widths of some resonances 
approach zero as h —> hb 0, [l3|, [l4| . Studies of specific examples show that when h = hb the structure has a stable 
wave guiding mode localized near the structure which cannot be excited by an incident radiation despite the fact 
that the spectral parameters of this mode lie in the spectral range of the scattered modes (the radiation continuum) . 
In other words, such wave guiding modes are a new type of localized solutions of Maxwell's equations which are 
analogous to quantum mechanical bounds states in the radiation continuum. 

It has been observed that in the spectral range near these localized solutions, there is a local amplification of near- 
fields (as compared to the amplitude of the incident radiation) for periodic scattering structures. The amplification 
effect becomes stronger if the system exhibits more narrow resonances 0, Il2l - fl4j . 

The similarities between resonant quantum and electromagnetic scattering phenomena are evident. It is therefore 
natural to develop a similar mathematical formalism of Siegert states to study the resonant electromagnetic scattering 
for periodic structures. Although a qualitative analogy between electromagnetic "trapped" modes and quantum 
Siegert states has been drawn, no rigorous quantitative studies have been carried out. It appears that there are two 
substantial differences between the theories of quantum mechanics and electromagnetism which prevent one from a 
trivial extension of quantum mechanical Siegert states to electromagnetism. 

First, the electric field E propagating in a non-homogeneous medium with a dielectric susceptibility e satisfies the 
wave equation -^d 2 E = AE. Consequently, the stationary analog of the Schrodinger equation, which is obtained 
when E(r,i) = e - iwt E w (r), is: 

-A + k 2 (l -e(r)) 

This equation can still be viewed as an eigenvalue problem, but in contrast to the stationary Schrodinger equation, 
the "potential" V(r) — fc 2 (l — s(r)) depends on the spectral parameter k 2 . Furthermore, the dielectric susceptibility 
could also be a complex-valued function of k if the medium is dispersive, i.e., e = e(k, r). 

Second, for a periodic scattering structure, the amplitude E w must obey Bloch's periodicity condition which is not 
compatible with the asymptotic boundary condition for quantum mechanical Siegert states. The stated differences 
require a modification of the very definition of Siegert states in electromagnetic theory. This problem is investigated 
in the present study. 

In Section UTJ electromagnetic Siegert states are defined for periodic scattering structures by means of the theory of 
compact integral operators when the polarization of the electromagnetic wave is preserved, i.e., when the vector wave 
equation is reduced to the scalar one. A relation between the Siegert states and the electromagnetic resonant scattering 
of waves is established. In Section ITVl the exponential decay of electromagnetic Siegert states excited by an incident 
wave packet is investigated. This allows one to identify the constructed Siegert states with " trapped" quasi-stationary 
electromagnetic modes observed in aforementioned numerical studies of periodic structures. In quantum mechanics, 
resonances of a system that consists of two coupled resonant systems depend on a coupling parameter h, and this 



E w = k 2 E^, k 2 = -j 
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dependence is determined by the theory of coupled resonances. A similar problem for electromagnetic resonances 
is studied in Section IHIl In the said section, the central result is an extension of the Kato-Rellich theorem [TtJ 
on analytic perturbation of resonance poles in quantum mechanics. In a similar fashion, the theorem established 
states analytic properties of electromagnetic Siegert states as functions of the coupling parameter h (the Regular 
Perturbation Theorem). By means of this theorem, it is then proved that if the scattering system admits a bound 
state in the radiation continuum (a resonance with the vanishing width) at some value of the coupling parameter 
h = hb, then the near-field can be amplified as much as desired by a suitable adjustment of the value of h near hb- 
This result provides a general mechanism to enhance optical nonlinear effects in periodic structures. Its efficiency 
has been recently demonstrated with a specific example to generate second harmonics by a double array of dielectric 
cylinders [l5[ (as depicted in Fig. QJb)). 



II. ELECTROMAGNETIC SIEGERT STATES 



Consider a scattering problem for an electromagnetic plane wave impinging a structure made of a non-dispersive 
dielectric. The structure is described by a dielectric function e(r) that defines the value of the dielectric constant at 
every point r occupied by the structure and it equals the dielectric constant of the surrounding medium otherwise. 
Without loss of generality, the surrounding medium is assumed to be the vacuum, i.e., e = 1. The material is said 
to be non-dispersive if its dielectric constant does not depend on the frequency of the incident wave. The structure 
is assumed to have a translational symmetry along a particular direction. In this case the dielectric function is 
independent of one of the spatial coordinates, say, the y coordinate. If the incident wave is polarized along the y— axis 
(the electric field is parallel to the y— axis), then the scattering problem can be formulated as the scalar Maxwell's 
equation: 

^d?E = AE, A = d 2 x + d* (1) 

where c is the speed of light in the surrounding medium (the vacuum), and E is the electric field. For a planar 
structure, the function e differs from 1 only within a strip (x, z) € (— oo, oo) x (—a, a). A planar structure is periodic 
if e(x + D g , z) = e(x, z) where D g is the period. In what follows, the units of length are chosen so that D g = 1. It 
is further assumed that the structure's dielectric constant exceeds that of the surrounding medium, i.e., e(r) > 1 for 
all r, and e is bounded. Two examples of such structures are shown in Fig. [T] (Panels (a) and (b)) depicting periodic 
arrays of infinite dielectric cylinders, all parallel to the y— axis. The function e(r) is piecewise continuous. For general 
periodic structures, e is invariant under the action of a suitable affine Coxeter group in the xy— plane. In this case, 
vector Maxwell's equations should be used as the polarization of the incident wave is not preserved in the scattering 
process. This general problem will not be studied here. 

If the incident wave has a fixed frequency, then the field E has a harmonic time dependence E(r,i) = E ul (r)e~ wt 
in Eq. (fTJ) where the amplitude E u satisfies the equation: 

A^(r) + fc 2 £ (r)^(r)=0, fc 2 = ^ (2) 

Let the x, y, and z axes be oriented by the unit vectors ei, e2, and e3, respectively. In the scattering theory, E u is 
sought as a superposition of an incident wave, and the corresponding scattered wave E^: 

E u (t) = e* k ' r + £^(r), k = fc x ei + fc z e 3 (3) 

where k is the wave vector of the incident wave. The units of the field E are chosen so that the amplitude of the 
incident wave is 1. The scattered wave obeys the outgoing wave boundary conditions at the spatial infinity. The 
periodicity of the scattering structure requires that the amplitude E^ satisfies Bloch's periodicity condition, 

E u {x + l,z) = e ik *E u (x,z) (4) 

Under the specified boundary conditions, the amplitude E u satisfies the Lippmann-Scwinger integral equation: 

E u (r)=tt[E u ](r)+e? k - T (5) 

where H is the integral operator defined by, 

H[£](r) = / (e(r ) - l)G fe (r|r )i?(ro)dro 
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FIG. 1. Panel (a): An example of scattering structures considered in this work. Infinitely long parallel cylinders are placed 
parallel to a y-axis periodically along an :r-axis in the vacuum eq — 1. The cylinders are characterized by a dielectric constant 
e c > 1. The rectangle D in Eq.© is enclosed by the dashed line, and the z-axis. 

Panel(b): Two periodic arrays such as the one on Panel (a) are placed parallel to each other at a distance 2h between the axes 
of any two opposing cylinders. 

Panel (c): The Cfc^ plane. The cuts run vertically from the diffraction thresholds (indicated by empty circles on the real axis) 
into the lower half of the complex plane. The diffraction thresholds divide the real axis into a countable set of intervals denoted 
I;, I > 0. The interval Iq lies below the radiation continuum. The rest of the intervals partition the radiation continuum. 



The kernel Gfc(r|ro) is the Green's function for the Hclmholtz equation, (A + k 2 )Gk(r\ro) = — 47T(5(r — ro), with the 
outgoing boundary condition. It is given by Gfc(r|ro) = inH$(k\v — y$\) where Hq is the Hankel function of the first 
kind of order 0. 

The existence of solutions to the Lippmann-Schwinger integral equation that also satisfy Bloch's condition is 
established by extending the operator valued function k 2 i— > H(/c 2 ) to a suitably cut complex plane, and thereafter, 
applying the Fredholm theory of compact operators. To this end, let the x~ component k x of the incident wave vector 
k be fixed. The magnitude k — |k| varies only in its z— component, k z = ±^k 2 — k 2 . Now let S E be the support of 
the function (x, z) — > (s(x, z) — 1) in the strip S = [0, 1] x (-co, oo) of the x, z-plane, and let D be a rectangle in S 
containing S e , i.e., 

S £ CD = [0,1] x [*_,«+] (6) 
Then for an amplitude E u satisfying Bloch's condition, 

B[E u ]{t) = ^Y, e " nkx I ( £ ( r °) ~ l)Gfc(r|r - me x )E u {v Q )dv Q (7) 

where 7L denotes the set of all integers. This equation defines H as an operator on L 2 (D), the Hilbert space of square 
intcgrable functions on D with respect to the Lebesgue measure. Note that H does not depend on the rectangle D 
which to some extent is arbitrary. Indeed, the integrals of Eq.(|7|) extend only over the support S E of the function 
[x, z) i — ^ [s{x, z) — 1) in the strip S — [0, 1] x (-co, oo) of the x, z-plane. The rectangle D is only introduced to obtain a 
connected and compact region of integration, which is convenient for the subsequent analysis. In general, the support 
of the function [x, z) t— > (e(x, z) — 1) is not necessarily connected as, for example, in the case of multiple scatterers in 
the strip S depicted in Fig. (2fb). The solution to the Lippmann-Schwinger integral equation ([5]) will be sought in the 
Hilbert space L 2 (D). Such a solution then extends naturally to the full x,z— plane by Bloch's periodicity condition 
(J4j> , and the Lippmann-Schwinger integral equation. 

It is not hard to see that the summation and integration can be interchanged in Eq. ([7]) even though the underlying 
series is only conditionally convergent. The Poisson summation formula is then applied to yield the following form 
for the integral operator H, 



H[£](r)= / (e(r )-l)H(k 2 :r-r )E(r )dr , E e L 2 (D) (8a) 
Jd 
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where the function H is given by, 



H(C;x,z) = —} , k x . m = k x + 27rm, C 7^ kl , m £ 7L (8b) 

I i —L r _ hi 



The square roots are defined by choosing the branch cut of the logarithm along the negative imaginary axis, i.e., if it) 
is a nonzero complex number, then 

111 ^77 , . 

log w = In \w\ + 1 argw, — — <argw<— (9) 

In the scattering theory, the branch points k 2 m , m £ %, are called the diffraction thresholds, and they will be referred 
to as such in what follows. 

Note that for w real, then Imfy 7 !/;} > 0. Therefore, Im | — k 2 m j > for all but a finite number of diffraction 

thresholds k 2 m such that ( = k 2 > k 2 m . It follows that in the series of Eq. (|8b[) all but a finite number of terms 
decay exponentially as \z\ — > 00. Hence the said series converges uniformly on D. The uniform convergence on D still 
holds when the range of the variable £ = k 2 is extended to the cut complex plane (Dfe x obtained by excluding all the 
vertical half lines running from the diffraction thresholds into the lower half of the complex plane, i.e., 

<C fe . E = <D - |J {kl m -is,s< 0} 

This is because in this case too, Im j ^Jc, — k 2 m j > except for a finite number of terms. It follows that £ H> H((; r) 

extends to an analytic function on (D^. Figure [TJc) shows a sketch of the cut plane C^. 

Since for all C € ( Cfe x , the kernel (r, ro) i-> (e(ro) — 1)H(C; r — ro) is square integrable on D x D, the integral operator 
H(C) is Hilbcrt-Schmidt in £{L 2 (D)} , the space of bounded linear operators on L 2 (D). In particular, H(C) is compact 
for each £ £ <Gk x ■ Moreover, the operator valued function defined on (Dfc^ by ( 1— )• H(^) is analytic. By the analytic 
Fredholm theorem [16], it then follows that if the inverse operator [1 — H(C)] -1 exists at some point ( £ <Ck x , it must 
be meromorphic throughout <Ck x ■ That the said inverse exists at some point is obvious. Indeed, as C - ^ in (Dfc x , the 
norm of the operator H(£), i.e., ||H(£)|| = ||(e(-) — 1)H((; •)IU 2 (Dx_d), also converges to zero. In particular for £ near 
0, ||H(C)|| < 1 and therefore [1 — H(^)] -1 exists as a Neumann series. Thus, the generalized resolvent £ H> [1 — H(£)] _1 
is meromorphic in <Dk m ■ 

The poles {Cn}n of [1 — H(C)] -1 are isolated and form a discrete set. The same analytic Fredholm theorem guarantees 
that at each of the poles £ n , there exists a nonzero solution to the generalized eigenvalue problem 

H(C« )[-£'«] = E n 

The generalized eigenfunctions E n will be referred to as the Siegert states. 

It will be assumed throughout the rest of the work that the Siegert states are nondegenerate, i.e., the poles ( n are 
all simple. The results to be derived in the rest of the work could easily be generalized to the case when the poles 
are not simple. Yet, to our knowledge, there has been no report of resonance poles of multiplicity greater than one 
in either the literature of electromagnetism or that of quantum mechanics. However, we do not have a rigorous proof 
that higher multiplicity poles cannot occur in the studied scattering problem. For the structures depicted in Fig. [I] 
perturbation theory (when the radius of cylinders is much less than the structures period) shows that all the poles 
are simple [l3|, 0] • 

The assumed nondegeneracy of the Siegert states {E n } n implies that the residues {H n }„ of [1 — H(£)] -1 at the 
poles {Cn}n are rank one operators on L 2 {D). In other words, if (/, g) = f D /(ro)g(i"o)dro is the inner product on the 
Hilbert space L 2 (D), then 

Vn, 3<p n £ L 2 (D) : W> £ L 2 (D), H„[V>] = (tp n ,ip)E n (10) 

Since the incident wave Ei(k 2 ;r) = e ikr is analytic in ( = k 2 on (Dfe^, it follows that the amplitude E u = [1 — 
K(()]~ 1 [Ei((; •)] extends to a meromorphic function of the variable £ = k 2 on (D^. Its partial fraction expansion 
reads 

E u (v) = e^ r + V 2 °" E n (r) + E a (k 2 ; r), a n = (<p n , E t (k 2 - ■)) (11) 

^— ' k z — kf +lln 1.2 _t.2 _»r 



k 2 =kl-iT„ 
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where each pole ( n has been decomposed in its real and imaginary parts as £ n = k 2 n — iT n owing to the fact that 
the imaginary parts of the poles are nonpositive as will be shown shortly. The remainder E a is analytic in k 2 on the 
cut plane C^. It describes the so called background or potential scattering, while the Siegert states account for the 
resonant scattering. 

Since the Siegert states are solutions to the homogeneous linear equation H(£)[_E] = E, they should be normalized 
in some way in order for the functions (p n in Eq. (|10j) and the coefficients a n in Eq.(|ll[) to be uniquely defined. A 
natural normalization condition will be introduced in Section 1111 Al for Siegert states that arise as perturbations of 
bound states in the radiation continuum. 

That all the poles of the generalized resolvent £ M> [1 — H(C)] -1 have nonpositive imaginary parts, i.e., T n > 0, Vn, 
follows from a series of observations. First, note that a Siegert state E n corresponding to the pole £„ = k 2 , — iT n 
satisfies Eq.J2]) for k 2 = £ n , and therefore 

~E~&E n - E n AK - 2tT n e\E n \ 2 = 

By Green's theorem, it follows that, 

2iT n [ |£ n (r )| 2 £(r )dr = / (£ tl (r )V£„(r ) - £„(r )V£„(r )) ■ dn (12) 

JD' JdD' 

where D' = [0,1] x [21,22] is any rectangle containing the rectangle D, and no is the outward normal. By Bloch's 
periodicity condition, the integrals on the line segments (x,z) g {0,1} x [21,2:2] of dD' are canceled out so that the 
integral on the right side is to be carried out only over the line segments (x, 2) g [0, 1] x {zi, 22} of dD' . The result 
can be expressed in terms of the scattering amplitudes associated with the asymptotic behavior of the Siegert 
state E n : 



\ x,m V sn x , m / z < Z- 



£ „ WH ™« z>zi (13.) 



Equations (|13a,|) follow immediately from equations ([8]). In particular, the amplitudes are given by the formulas, 

S± = , = — / (e(r ) - l)g tt (r K ( -*° fc ^W^- fc ^)rfr , r = x ei + 2 e 3 (13b) 
2./Cn - k 2 Jd 



(14) 



Evaluating the right hand side of Eq. (fT2")) yields, 

_ E mez Rc (\s+fe-^iV^^} + \s m \ 2 e 2 ^{V^^} 

J D , |i;„(ro)| 2 £(ro)dro 

The series in Eq. (|14p is then split into the sums over the two complementary index sets, 

r*(Cn) = {m € TL I Re{C„} > k^ m }, I ci (C„) = {m g % \ Re{C„} < k^ m } (15a) 

In the scattering theory, when £ = k 2 is real, the sets I op (fc 2 ) and l cl (k 2 ) label open and closed diffraction channels 
respectively; thus the superscripts "op" for open, and "cZ" for closed. 

The choice of the branch cut for the logarithm given in Eq. © ensures that 



Im 



yCn~kl m } >0, VmeI d (Cn) (15b) 



In particular, in the series of Eq. (|14[) . the contributions from terms whose indices lie in I c '(£ n ) decay exponentially as 
22 — > +00 and 21 —¥ —00. In the said limits, the aforementioned series is therefore reduced to a finite sum over the 
set Pf(C„): 

E m ei°HO Re { ~ kj m } (\S+\ 2 e- 2 ^{V^^} + \S m \ 2 e 2 ^{V^^}) 
21^-^ J D , \E n (r )\ 2 e(ro)dr 

Observing that Re | sfCn ~ k 2 m | > for all m g I of> (C„), and since the dielectric function e is positive, it follows 

that r„ > 0, and therefore the imaginary part of the pole is indeed necessarily negative or zero, i.e., all the poles 
are in the lower half of the cut complex plane (Dfe x . 
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III. REGULAR PERTURBATION THEORY 

Suppose that the dielectric function e depends on a real parameter h which is called a coupling parameter. The 
simplest example of such a coupling is given by two periodic arrays of dielectric scatterers that are parallel and 
separated by the distance 2h (Fig.QJb)). The arrays are embedded in a medium of dielectric susceptibility 1. Each 
array is characterized by a dielectric function Ei(x,z), i = 1,2 such that e$ > 1 on the scatterers, and £j is of period 
1 in the a;-direction. The resulting dielectric function describing the scattering of light on the structure is then, 

e(h;x,z) = 1 + {e\(x,z — h) — 1) + (e2(x, z + h) — 1) (16) 

In the most general case, the dependence of the dielectric function e on the coupling parameter h implies that the 
poles of the generalized resolvent £ h->- [1 — H(/i,£)] _1 also depend on h. If the coupling is sufficiently smooth as 
indicated in the theorem stated below, then the poles of the generalized resolvent as well as the corresponding Siegert 
states depend continuously on h. 

Before stating the theorem, a clarification must be made on the boundary of the rectangle D in Eq.©. This 
rectangle was chosen to contain the support S E of the function (x, z) h-> (e(x, z) — 1) in the strip S = [0, 1] x (— oo, oo) 
of the x, z-plane. In the current situation, the dielectric function depends on the coupling parameter h, and therefore 
the support Sw/,) in question could change with h resulting in a different choice for the set D. This is the case, for 
instance, in the example of the two parallel arrays separated by the distance 2h. If the two scatterers in the strip 
S are taken further apart, the rectangle D is stretched further accordingly. So, in what follows it will always be 
assumed that h varies in an open possibly finite interval Jo such that U/igj ^e(fc) is bounded, and the rectangle D 
will be chosen to contain the latter set. The Theorem on the regular perturbation of electromagnetic Siegert states is 
formulated as follows: 

Regular Perturbation Theorem. Suppose that the map (h, £) i-> H(/i, £) from JoxC/^ to £{L 2 (D)} is continuously 
Frechet differentiable. Further, suppose that for some h n £ Jo, the generalized resolvent ( i-4 [1 - H(/i n ,£)] has 
a simple pole £ <Ck x , and let E n be the corresponding Siegert state. Then there exists an open interval J C Jo 
containing h n , and a unique continuously differentiable function h i— > £ n (h), h £ J, such that Cn(h n ) — Cn, and 
for all h £ J, (n(h) is a simple pole of the generalized resolvent £ — > [1 — H(ft, C)] _1 - If E n (h) is the Siegert state 
corresponding to the pole ( n (h), then the function h > E n (h) from J to L 2 (D) is continuously differentiable. 

This theorem is an extension of the Kato-Rellich theorem [l7| to the present case. Due to its length and complexity, 
the proof is left to the Appendix. It should be noted, however, that if the functions (x, z) £i(x, z), i = 1, 2 in Eq. p^)) 
are piecewisc differentiable, then (h, £) i-> H(/i, £) is Frechet differentiable so that the theorem does indeed hold for 
two parallel arrays separated by the distance h. 



A. Perturbation of bound states in the radiation continuum 

The diffraction thresholds k 2 . m on the real line depend on the x— component k x of the incident wave vector. They 
can be ordered independently of the parameter k x by defining a sequence Q_ m {k x ) = (27rm± |[fc x ]|) 2 , m £ {0, 1,2, . . .}, 
where [k x ] is the argument of e lfex in the interval (— ir, ir]. The sequence {C± m (kx)}m=Q coincides with the sequence 
of diffraction thresholds {(k x + 2itm) 2 } ml z%, and for all k x , 

Co(M < CM < Ci(k x ) < CM < <2(kx) < C 3 (kx) < ... 

The threshold Cofc) i s called the radiation continuum threshold, and the interval Io = (— oo,Cq(^x)) is said to lie 
below the radiation continuum. This is because whenever k 2 < Co(k x ), then the scattered amplitude in Eq.Q 
necessarily decays exponentially in the spatial infinity, \z\ — > oo, as can be inferred from Eqs.Q, ([5]), and ©. Hence, 
no electromagnetic flux is carried to the spatial infinity in this spectral range. In contrast, if k 2 > £q (k x ), the amplitude 
E^ oscillates and represents outgoing (scattered) radiation that carries an electromagnetic flux to the spatial infinity. 
The spectral range above Co(^x) on a real line is therefore referred to as the radiation continuum. It is the disjoint 
union of the intervals, 

ii = (£(**), Ci (*«)), i2 = (Ci(**), £(*»)), i3 = (cr(M,£ 2 (M), ■•■ 

Note that when [k x ] is or tt, some of the intervals I; are empty, owing to the fact that some of the diffraction 
thresholds fuse. 
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Returning to Siegert states, recall that these states are generalized eigenfunctions to the generalized eigenvalue 
problem H(£„ )[E n ] = E„ on L 2 (D) with complex eigenvalues £„. These states are naturally extended to the whole 
x, z— plane by means of Eqs. (fT3|) . In general, if the generalized eigenvalue Cn — k„ — iT n lies below or in the interval 
Io in the (D/^ -plane, then the corresponding Siegert state is square integrable on the strip S — [0, 1] x (—00,00) of 
the x, z-plane. This is because for such states, the set I op (Cn) of Eq. (115a[) is empty, and therefore by Eqs. (|13a[ ) and 
(|15b|) . these states decay exponentially in the asymptotic region \z\ — > 00. In particular, the Siegert states for which 
the pole £{, is real and less than the continuum threshold Co(^x) are the bound states below the radiation continuum 
of the system. 

On the contrary, the Siegert states E n whose corresponding generalized eigenvalues Cn — fe„ — *r„ lie below an 
interval Ij, i > 1 of the radiation continuum, i.e., k 2 € I; and T n > 0, are not necessarily square integrable on the strip 
S. This is because for these states, the set I op (( n ) is not empty, and consequently, the terms of the series in Eq. (|13a[) 
indexed by this set are unbounded. However, when the pole £ n is in Ij, J > 1, i.e., T„ = 0, it will be shown shortly 
that the resulting Siegert state is square integrable on the strip S and therefore, such a Siegert state is a bound state 
in the radiation continuum. The rest of this subsection is devoted to the proof of this assertion. 

To proceed, suppose that for some value hb of the coupling constant h, there exists a Siegert state E\, whose 
corresponding generalized eigenvalue Cb is real and lies in an interval I;, I > 1, of the radiation continuum. By 
the Regular Perturbation Theorem there exist continuously differentiable functions h i-> ( n (h) and h H> E n (h) on an 
interval J containing hb such that C, n {hb) = Cb, and E n (hb) = Eb- Furthermore, E n {h) is the Siegert state corresponding 
to the pole Cn(h>) of the generalized resolvent C H> [1 — H(/i, C)] -1 for all h G J. Put ( n (h) — k 2 (h) — iT n {h) as before. 
Then T n (hb) = 0. Without loss of generality, the interval J is assumed to be sufficiently small so that for all 
h E J\{hb}, r n (/i) 7^ 0. Equation (TlH) for the Siegert state E n {h) is then rewritten as 

/ \E n (r )\ 2 e(h;r )dr = -L £ Ref/^u} (|S+| V 2 WV£^} + \S^e^{V^^}) 

where it is understood that the Siegert state E n , the pole Cm and its imaginary part — iT n , as well as the amplitudes 
Sm are all functions of h. Note that the values of z\ and z-i are assumed to be independent of the coupling parameter 
h. This is because h varies in a small interval J, and therefore the values of z± in Eq.Q, which now depend on h are 
bounded. The condition z\ < z~(h) < z + {h) < z%, V/i € J, can therefore be realized by choosing z\ and Z2 sufficiently 
large. 

As h — > hb, then E n — > Eb in L 2 (D), and therefore the left hand side of Eq. lfTTj) remains finite as the dielectric 
function e(h; •) is bounded. The limit of the right hand side as h — > hb may also be calculated by first computing 
explicitly the complex square roots involved according to the logarithmic branch cut described in Eq.(|9|). Put C~m{h) = 
k 2 (h) -k^ m ,m€%. Then 



\/(n(h) - k% 



if me l op ((n(h)) 



2 ( y/u (h) 2 +r„ (h)* -u (h)) 



where l op and l cl are the index sets of Eq. (|15ap . In particular, if m £ I op (£ n (/i)), then £, m {h) > 0, whereas £ m (h) < 
whenever m £ I cl (( n (h)). As h — > hb, T n (h) —> 0, and, 

/ \E b (r )\ 2 e(hb;r a )dr = lim — L- £ y/tJKj (\S+ (h)\ 2 + \S m \ 2 

In particular, for each m € I op (£b), the limits lim/j-^,, ^" Wr must exist and must be finite. By continuity of the 
functions h H> Si^(h) on J, there exist complex numbers b such that 



^ "tI^t = S ^ Vm e rP(Cfc) (18) 
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As z\ — ¥ — oo and Z2 — > oo, the second series in Eg. (IT7j) converges to zero. Therefore the function r — > E},(r)\J e{hb\ r) 
is square integrable on the strip S = [0, 1] x (-co, oo) of the x, z— plane, and 

f \E b (r )\ 2 s(h b ;r Q )dr = £ 6 | 2 + 1^/) (19a) 

Js ™ei°p(c„) 

Since, e(h b , ■) > 1, it follows that I?;, e L 2 (S) as claimed. Thus, if is a bound state in the radiation continuum, it 
can always be normalized by the condition 

/ \E b (r )\ 2 e{h b -rv)dr = 1 (19b) 
Js 

This result justifies the term "bound state" introduced by analogy with quantum mechanics where bound states 
represent square integrable eigenfunctions of the Hamiltonian operator, i.e., an electromagnetic bound state is a 
localized solution of Maxwell's equations. 

Finally, by continuity E n (h) — >• E b as h — > h b , and therefore the normalization condition (|19bj) determines uniquely 
the Siegert states E n {h) along the curve h*-> ( n (h), he J. 



B. Near field amplification mechanism 



Resonant phenomena in the scattered electromagnetic flux can be described by the formalism of Siegert states as 
given in Eq. (lll[) . However, a complete description requires calculating the residues a n in Eq. tjlip . Here these residues 
are calculated for Siegert states that arise as perturbations of bound states in the radiation continuum in the sense 
of Section Till A\ i.e., for the states E„ (h) where \h — h b \/h b <C 1. In particular, it is shown that if the scattering 
structure has bound states in the radiation continuum, then there exist regions (the so called "hot spots") in which 
the field amplitude can be amplified as much as desired by taking the value of the coupling parameter h close enough 
to h b . The unbounded local growth of the field amplitude is essentially due to the linearity of Maxwell's equations. It 
disappears when a non-linear dielectric susceptibility, required for large field amplitudes, is included into Maxwell's 
equations [l5j]. Nevertheless, such a local field amplification enhances quite substantially optical nonlinear effects in 
the scattering structure as shown in 15|. The following analysis also suggests that the concept to enhance optically 
non-linear effects by using bound states in the radiation continuum is universal because the existence of "hot spots" 
is proved to be a characteristic feature of such scattering structures. 

When estimating the coefficients a n , it is convenient first to give another equation for them that is alternative to 
that of Eq. dlip . Next, this equation will be analyzed in the vicinity of a bound state in the radiation continuum. In 
particular, if a bound state in the radiation continuum E b exists at the critical value h b of the coupling parameter, 
and J is the interval produced by the Regular Perturbation Theorem, while E n (h;r) is the Siegert state that arises 
as a continuous perturbation of the bound state E b for h 6 J, then it will be proved that a n (h) ~ y/T n (h) as h — > h b . 
Recall that £ n (h) = k 2 (h) — iT n (h) is the generalized eigenvalue at which the Siegert state E n (h) exists. 

To proceed, put E u = E u where E u is the amplitude in Eq.([3]). Then by the decomposition (fTTj) . E u — > E n as 
C — > Cn< From the system, 

AiL + Ce{h)E u = 
AE n + (e(h)E n = 



it is derived by Green's theorem that, 



/ (E n VE u - E u VE n ) ■ dn + / E u E n e{h)dr = 

JdD' V ' JD' 



where D' is the same rectangle as in Eq. (| 12|) . Then, by splitting the amplitude E u in terms of the incident and 
scattered waves, E u (r) = ~ e tk ' r + -E 1 * (r), one infers that 

_ (C - Cn) U aD ,(ikEZ - VK)e* k - r ° ■ dn + (C - C») J D , E~e^e(h)dr ] 



J dD , (EnVE' - E'VE n ) ■ dn + (C - C«) f D , E^E n e(h)dr 
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where it is understood that E n and £„ 
is obtained by taking the limit £ — > Cn 

a n (h) 



are functions of h. Since a n is independent of £, the desired expression for a r< 



J dD , (ik n E n - VE n ) e ik »' r ° • dn 



low 



E n Vd c E*, - d c EtVE n ) ■ dn 



X 1 



|£ n | 2 e(^r 



k n — k 



C=Cr, 



(20) 



where for the first term in the denominator, I'Hopital's rule has been applied. This formula would not generally be 
useful as the term E^ involves the coefficients a n implicitly. However, if the Siegert state E n (h) is taken near a bound 
state in the radiation continuum, as assumed here, the first-order perturbative expression of a n only depends on the 
Siegert state E n (h). 

Such a perturbative expression is obtained by analyzing each of the integrals in Eq. (l20l) separately. First, it is 
observed that as h — > hb, then E n {h) — > Eh- By letting z\ — > — oo, and z 2 — > oo, it follows from the normalization 
of Eq. (|19bp that the second integral in the denominator of Eq. (|2"0)) can be made arbitrarily close to 1 provided h is 
sufficiently close to hb- Hence, if the first integral of the said denominator can be shown to converge to zero in the 
limits considered, this would imply that in the leading order of perturbation theory, the denominator of Eq. (|20|) is 
approximated by 1. That the first integral of the said denominator does indeed converge to zero as h — ¥ hb, z\ — > — oo, 
and Z2 — > oo can be established through a lengthy, but straightforward calculation whose details are omitted. It 
can be carried out along the following lines. As noted earlier, Bloch's condition implies that the first integral in the 
denominator of Eq. (|2"0|) is to be carried out on the segments {x, z) € [0, 1] x z 2 } of the boundary of D' . Now, 
as h — > hb, the Siegert state E n (h) becomes a bound state Eb, and therefore, it decays exponentially in the spatial 
infinity \z\ — > oo. Similarly, the same exponential decay can be established for the terms involving the derivatives of 
E^ in the limit h —¥ hb- It then follows that in the limits h — > hb, Z\ — > — oo and z 2 — > oo, the integral in question 
vanishes. Thus the denominator of Eq. (|20| remains indeed close to 1, provided h is sufficiently close to hb. 

On the other hand, the numerator to Eq. (f2"U| can be evaluated in terms of the amplitudes in Eqs. ([T3|) to yield 



/ (ik n E n - VE n ) e lk "' ro • dn = ^/T n (h)A n (h, Zl ,z 2 ) 
Jd 1 



where if k n = k x B\ + yCn — k^e^, then 

1 J v r «(") */r„(/).i 



The case k„ = k x ei — y/C, n — fc^ e 3 is similar. As h — > hb, then Cn{h) — > (b, and since > k x , it follows that 
^y\/Cn{h) — fefj — » 0. Hence, in the said limit, A n (h, zi,z 2 ) — > 2iS^ b ^/C,b — k x for S^ b given in Eq. flTBl) . Thus, in 



Im 



general, a n can be written as 



a n {h) = a n {h)y/T n {h) 



(21) 



where a n (h) is bounded and has the property a n (h) — > 2iSQ b ^/(b — k x as h — > hb- 

The most important consequence of the structure of the coefficients a n near a bound state in the radiation continuum 
is a local amplification of the amplitude E u as compared to the amplitude of the incident radiation. Indeed, if 
the wavenumber k of the incident radiation and the coupling parameter h are tuned to satisfy the condition k = 
k n (h), V/i G J, then the amplitude (TTT1) becomes: 



where h i— > \\E a (h, -)||l 2 (_d) is bounded on J. It follows that, 

\a n (h)\ 



(22) 



> 



y/T^h) 



E, 



n\\L2(D) 



\E 



a\\L 2 (D) 



As h — > h b , then E n — > E b in L 2 (D), and therefore ||£7 n ||i, 2 (r>) — ► II^Hl^d) 7^ 0- Hence if the constant Sg b in Eq. (j2"Tj) 



is nonzero, then, 



lim \\E u \\ L 2r D) = oo 

h—thb 
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and this can only happen if the amplitude E w diverges in some regions of the rectangle D. Even though the fact that 
a n (h) is always nonzero in the limit h — > h b (i.e., Sq b 7^ 0) could not be verified for all kinds of couplings, there are 
many systems in which it holds true. For example, this is the case of the normal incidence (k x = 0) for a symmetric 
double array depicted in Fig. [TJb) when the dielectric functions e\ and £2 in Eq. flTrH) are identical, and symmetric 
with respect to the reflection (x, z) H> (x, —z). In this particular case, the parity operator, S[E](x, z) = E(x, —z) on 
L 2 {D), and the Lippmann-Schwinger integral operator H commute so that Siegert states are always symmetric or 
skew symmetric with respect to the reflection (x, z) 1— > {x, —z) in the x, z-plane. In particular, the amplitudes 
given in Eqs. (fT5)) are such that = ±5~ depending on whether the Siegert state they correspond to is symmetric 
or skew symmetric. It follows that if the bound state E b happens to lie in the interval Ii of the (Dfc x -plane (case of 
one open diffraction channel), then Eqs. (TT9|) are reduced to 



2va 



so that the coefficient in question is indeed nonzero and the amplitude E u is amplified in the vicinity of the bound 
state in the radiation continuum Eh- Note that the amplification of the amplitude E^ was also observed perturbatively 
in the case of a more general incidence angle (i.e., k x 7^ 0) when the bound state E b lies in the intervals Ii and I2 of 



the radiation continuum 14 1 



The aforementioned amplification can only happen in a region near or within the scattering structure. This can be 
understood by analyzing the first term of Eq. (|22[) from which the amplification should result. Outside the scattering 
region, the Siegert state is expressed in terms of its scattered amplitudes by Eq. (|13a[) . As before, the series involved in 
this expressions split over the index sets I op (( b ) and I cl {( b ) ( for h near h b , I op (C„(/i)) = I op (C&) and I cl (( n (h)) = I cZ (C&))- 
On one hand, the terms of the series whose indices lie in I c '(Cf>) decay exponentially in the spatial infinity so that 
no amplification of the field can be obtained from them. On the other hand, the terms whose indices lie in I op ((f,), 
disappear near the bound state in the radiation continuum as these terms are proportional to y/T n (h) as indicated by 
Eq.(|18p. From the physical point of view, the noted local field amplification results from a constructive interference 
of scattered fields from each scatterer (an elementary cell) of the periodic structure. It is important to note that the 
amplification magnitude can be regulated by varying the coupling parameter h, which provides a natural physical 
mechanism to control optical nonlinear phenomena if the structure has a nonlinear dielectric susceptibility. This 
mechanism has been used to demonstrate that a double periodic array of dielectric cylinders (depicted in Fig. [T](b)) 
with a non-zero second-order nonlinear susceptibility can convert as much as 44% of the incident flux into the second 
harmonics when the distance between the arrays is properly adjusted. The smallest distance at which such a high 
conversion rate can be achieved is about a half of the wave length of the incident radiation [l5[ . 

IV. DECAY OF ELECTROMAGNETIC SIEGERT STATES 

A Siegert state can be excited by an incident radiation, e.g., by a wave packet whose spectral distribution peaks 
around a desired wave number k c w k n . When passing through the structure, some of the wave packet energy is 
trapped by the structure and remains in there for a long period of time, decaying slowly by emitting a monochromatic 
radiation. This is established in a fashion similar to that of the study of the decay of unstable states in quantum 
mechanics [l|. To illustrate the principle, only the case of normal incidence is considered here (i.e., k x =0). Other 
incidence angles can be treated in a similar manner. 

The exact statement which will be proved is as follows. Suppose that a Siegert state E n exists at the pole £ n = 
fc„ — iT„ such that k n > and r n > (i.e., the Siegert state is not a bound state in the radiation continuum). Then 
this state will decay with time by emitting a monochromatic radiation at a wavenumber k n , and with a life time r n : 



i n = M±Vk[+n — — (23) 

V 2 r„<si cT n r„<i cT n 



As a starting point, consider an incident wave packet, 

Ei(r,t) = / A(k) cos(k • r — u)t)dk, u> = ck 
Jo 



where A(k) is the distribution of wavenumbers in the wave packet. Then the solution E of the wave equation (p} 
reads 



S(r,t) = Re|jT A{k)E L0 {v)e- luJt dk^ 
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where E u is the amplitude of Eq. ([2]). From the meromorphic expansion of E u in Eq. (fTTj) it then follows that, 

E{v,t) = Re !^Ei(r,t) + ^ E n (r)Q n {t) + jf A(fc)^ a (fc 2 , r) e - llJJt dfc| 
where the time dependence O n (i) of the Siegert state E n is, 

n n (t) = s„Vr^ /°° r 5 4ri t e " cfc * dfc ( 24 ) 

JO _ ^ + !l „ 

for a„ defined in Eq. (|2ip . As t — > oo, then Q n (i) — >• as required by the Riemann-Lebesgue lemma. For a typical 
physical wave packet, the function A(k) decays fast as k — > oo and is also analytic in the complex k— plane. These 
properties allows for evaluating the integral (|24l) by the standard means of the complex analysis. It is also worth 
noting that bound states in the radiation continuum cannot be excited by an incident radiation because T n = 0. 

To avoid excessive technicalities of the general case, a specific form of the distribution A(k) is chosen to illustrate 
the procedure, which is sufficient to establish the main properties of the decay of Siegert states. The simplest form 
that is also most commonly used in physics is a Gaussian wave packet centered around a wavenumber k c : 

A(k) 



aV2w 

In the limit a —> 0, A(k) — > S(k — fc c ), and the monochromatic case is recovered. The analysis will be carried out for 
a Siegert state E n corresponding to a pole £„ = k\ — iT n such that A; 2 , > and r„ > 0. Siegert states with fc 2 < 
are discussed at the end of this section. 

The integrand in (|24[) is highly oscillatory for large t, the contour of integration must be deformed to a curve along 
which the phase is constant to obtain a fast converging integral, and thereby to determine the leading term. To this 
end, first, the change of variable £ = fc ~kc+?_ ctcr is made to obtain the following more amenable form for fl n : 

= o 2 7= / Tc — T^Tc — 1T\ ?± = 7k ( 25 

where the contour of integration Cq is a horizontal ray outgoing from the point £o = ^ k "^^ a toward the infinity in 
the £-plane as shown in Fig. [5J The position of the poles £± indicated on the same figure follows from the fact that 
for fc 2 > 0, then Re{v^} > while Im{VCT} < 0. 

Put (u,v) = (Re£,Im£). Since the function £ i-> decays exponentially in the region Re{£ 2 } > 0, the contour 
of integration in Eq. ([23|) can be deformed to a contour that consists of the constant phase curve C% = {(u,v) | uv = 
— u < Re£o} extending from £ to — oo and another constant phase contour C 2 : I m {£} = from — oo to oo. 
Figure [5] shows the modified contour. By Cauchy's theorem, VL n becomes the sum of three terms: 

n n (t) = a n (A n (t) + B n (t) + C„(t)) 

where A n is the residual contribution at the pole £+, while B n and C n are the contributions of the line integrals 
along the contours C\ and C2, respectively. It is proved shortly that it is the residual term A n that accounts for an 
exponential time decay of the Siegert state E n . For a long period of time, this term dominates in the decay radiation 
of a Siegert state, and it is only after it has decayed considerably that the term B n becomes dominant. The term C n 
remains small in comparison to A n and B n . 

To begin, the residual term A n is evaluated to yield the formula, 

A n {t) = m ^A(^)e- lct ^ 
V vi 

In particular, the wavenumber k n , and the life time r„ of the Siegert state given in Eq. (l23[) follow from the formula 
f° r VCn which is, 



\Jk\ — iT ri 



+ v^n + r 2 . r„ 



^2 (ki + ^WT^ n 
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Note that for T n <C 1 (a narrow resonance), the amplitude of the Siegert state is proportional to A(y^) A(k n ). 
Hence, the Siegert state corresponding to the scattering resonance at k = k n is excited if the Gaussian wave packet 
has a sufficient amplitude at k = k n , or, ideally, is centered at k n to achieve the maximal effect. 

The terms B n and C n also decay necessarily in time. For C n , the integral along C2 will necessarily remain bounded 
in time so that the multiplicative factor in Eq.(|25j) is dominant, i.e., 



C n (t)=0(e 



The same relation holds for B n , provided that 



) 



ct < 



(26) 



When this condition is violated, the point £0 moves into the region Re{£ 2 } < 0, and the exponential e - ^ starts to 
grow large on a part of the contour C\. Nonetheless, the term B n {t) still decays to zero in time. This is because, 
as noted earlier, the Riemann-Lebesgue lemma ensures the decay of O n (i) in time, and the terms A n and C n have 
been proved to converge to zero. For the actual decay rate of B n , it can be inferred from Laplace's method for the 
asymptotic expansion of integrals that B n = 0{l/t a ) for some positive number a, which is of no importance here 
however. Indeed, it is obvious that as t increases, the term B n should be larger than the residue at £ + , since the arc 
of the curve C\ in the region Re{£ 2 } < grows longer, and therefore the exponential decay of the Siegert state can 
no longer be observed above the background described by B n . The condition (|2"B"|) shows that the smaller the width a 
of the wave packet is, the longer the exponential decay can be observed. For instance, one can ensure the observation 
of the decay radiation of the state E n by requiring that the half-life time occurs before the decay amplitude becomes 
smaller than that of the background, i.e., 



cr n < 



h^ki + ^kjffl < r 



In other words, the width a of the packet must be comparable to VT^, if not smaller. 




FIG. 2. The £-plane for the integral in Eq. (|25|) . The region Re{£ 2 } > is shadowed. The exponential is bounded in 
it. The sketch is realized under the condition (|26[1 . when the point £0 is still in the shadowed region. When the condition is 
violated, part of the curve Ci lies in the region Re{£ 2 } < 0. 



For Siegert states with k 2 , < 0, the pole £ + is no longer enclosed by the contours Co, C\, and C2 if t is sufficiently 
large. The pole moves to the left of the contour C\ and the exponential time decay is absent, owing to the fact that 
the incident wave packet does not contain any radiation at wavenumbers close to k n for the Siegert state to be excited. 

Appendix: Proof of the Regular Perturbation Theorem 

Recall that V(/i, C) € Jo x € kx , H(/i, C) € £{L 2 (D)}, and that the map 

Jo x € kx -> £{L 2 (D)} 



(ft,0->H(ft,0 
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is continuously Frechet differentiable. For (ft, £) G Jo x <D& a , let R\(H.(h, Q) = [A — H(ft,C)] \ A G (D, be the resolvent 
of H(ft, C) when it exists, and let /o(H(ft,C)) = {A E (D|-Ra(H(/i, 0) exists} be the resolvent set of H(ft, C). The proof 
of the Regular Perturbation Theorem is aided with the following lemma: 

Lemma. Let S = {(ft, C, A) G Jq x (D^ x CD | A G p(H(ft, ())}■ The set S is open in Jo x (D/^ x (D cmc? nonempty. 

Proof. That 5 is nonempty follows from the results of Section [TTJ Indeed, if C is not a pole of C i-> [1 — H(ft, C)] _1 , 
then (ft, C) 1) G 5. To show that S is open, let (fto,Co, Ao) G 5. Then one has to prove that there always exists a 
neighborhood of (fto, Co, Ao) in Jo x <D/c x x <D that lies in S. 

Since p(H(ho, Co)) is open in <D, there exists Ai G /o(H(fto, Co)), Ai 7^ Ao so that [Ai — H(fto,Co)] _1 exists. Let then 
^ be the function, 

* : Jo x <£ k ^C{L 2 (D)} 
(ft,C)^ Ai-H(/i,C) 

The map VP is continuous, and ^(/lo, Co) is invertible. As the set of all invertible operators in £{L 2 (D)} is open, there 
exists an open neighborhood U C Jo x (Dk x of (ho, Co) such that V(ft, C) G U, ^(h, C) is invertible. 
Now by the first resolvent formula, 

i? Ao (H(/ lo ,Co))-i?A 1 (H(/ lo ,Co)) = (Ai -Ao)i?A (H(^o,Co))i?A 1 (H(^ ,Co)) 

so that, 

i? Ao (H(ft , Co)) (l - (Ai - X )Rx 1 (H(fto, Co))) - Rx, (H(ft , Co)) 

It follows that 1 — (Ai — Xo)R\ 1 (H(fto, Co)) is invertible. Now let $ be the function, 

$ : U x € -)• £{L 2 (L>)} 
((/ l ,C),A)^l-(A 1 -A)i? Al (H(ft,C)) 

Then $ is continuous. Since $(/io,Co, Ao) is invertible, and again, the set of all invertible operators in £{L 2 (D)} is 
open, it follows that there exists an open set W C U x (D such that for all (ft, C, A) G W, the operator &(h, (, A) is 
invertible. 

Thus, for (ft, C, A) G W , both Ai — H(ft, C) and 1 — (Ai — \)R\ 1 (R(h, C)) are invertible. Their composition is therefore 
invertible, and this composition is, 

(Ai - H(ft, C))(l - (Ai - A)i? Al (H(ft, C))) = A - H(ft, C) 
It follows that W C S. As IF is a neighborhood of (fto, Co, Ao) in Jo x (D/^ x (D, the set S is open. □ 
The proof of the Regular Perturbation Theorem is as follows. 

Proof. Let Cn be a simple pole of C ^ [1 — H(h n , C)] _1 for some fixed h n G Jo, and let E n be the corresponding Siegert 
state. Then 1 is an eigenvalue of H(ft,„,Cn) corresponding to the eigcnfunction E n in the usual sense. The theorem 
amounts to proving the existence of a curve h <— > C(^) m Jo x (Dfc^ along which the family of operators H(ft, C(^)) still 
has 1 as an eigenvalue. The corresponding eigenfunctions will then be the Siegert states of the said operators along 
the curve in question. 

The proof starts by providing a general formula for the eigenvalue Ao(ft, C) °f H(ft, C) which is the perturbed value 
of the eigenvalue 1 when the point (h n , C„) is displaced to (ft, C) in the Jo x (Dfc x space. Then by the Implicit Function 
Theorem, a curve h i-> C(ft) is found along which the eigenvalue X(h, C) remains 1, i.e., X(h,C(h)) = 1. 

As a starting point, note that as H(h n ,( n ) is a compact operator on the Hilbert space L 2 (D), the Riesz-Schauder 
theorem [l6[ implies that 1 is necessarily an isolated eigenvalue. Therefore, there exists 5 > such that 1 is the only 
eigenvalue of H(ft„,C n ) in the disk {A G (D | |A — 1| < 5}. It follows that V 5 = {(ft»,Cn,A)|A G (D, |A - 1| = 5} C S} 
where S is the set of the previous lemma. Since S is open, and V$ is compact, there exists an open W in J x (D/^ x (D 
such that V$ C W C S. Therefore there exists a connected neighborhood U of (h n ,( n ) in Jo x (D/^ such that 
A G p(R(h, 0) for all (ft, C) G U, and for all A G C with |A - 1| = S. 
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Put 

P(/*,0=o^/ [A-H^C)]" 1 ^, (h,()eU (27) 

Then P(/i, C) € £{L 2 (£>)}, and by Theorem XII.6 of [13, P(/i, C) is a projection. As the map (h, C) >-> P(h, C) of J7 to 
£{L 2 (Z?)} is continuous, it follows that the dimension of the range of P(h, C) is constant throughout Z7. 

Since the eigenvalue 1 of H(/i„,£„) is nondegenerate, it follows that the range of P(/i„,Cn) has dimension 1, and 
therefore the dimension of the range of P(h,0) is 1 throughout U. By Theorem XII.6 of [13, it follows that for all 
(h, C) G f, there exists a unique eigenvalue Ao(/i, C) 01 H(ft, C) m {A G (D | |A — 1| < 5}, and P(/i, £) is the projection 
on the corresponding eigenspace. In particular, P(/i, C)[-EVi] is in the eigenspace of Ao(/i, C), and therefore, 

K(h,()[P(h,C)[En}] =A (/i ) C)P>,C)[£n], (ft, C) 6 17 

Hence, 

C) fp(/»,C)[^]l) 

A (/i, C) = % O^u 

(E n ,P(h,0[E n }) 

where once again, (•, •) is the inner product on L 2 (D). This is the formula for the perturbed eigenvalue announced in 
the preamble of this proof. 

Now, observe that since (h, C) ^ H(/i, C) is continuously Frechet differentiable in U, so is (h, Q >-t P(h, (). Therefore, 
Ao is continuously Frechet differentiable in U. Now, Ao(/i ra , Cn) — 1, and as it will be shown shortly, d^Xo(h n , Cn) 7^ 0. 
It follows by the Implicit Function Theorem for Banach spaces that there exists an open interval J C Jo containing 
/i„, and a unique continuously differentiable function h ( n {h), h G J, such that (, n (h n ) = ( n and A (/i, CnC 1 )) = 1 
for all h £ J. Thus the proof of the Regular perturbation theorem is complete, provided the property d^Xo(h n , Cn) 7^ 
is established. The latter is achieved by studying the analyticity of the resolvent R\(K(h n , £)) = [A — H(/i„, C)] -1 in 
the variables A and £ separately when h = h n is fixed. 

Let A G C and C G Cfe* such that |A - 1| = S, (h n ,Q G U, and C 7^ Cn- Then the resolvent R\(R(h n , £)) = 
[X — H(h n , C)]^ 1 exists by definition of the neighborhood [/. Also, the resolvent Ri(H(h n , £)) = [1 — H(h n , C)] -1 exists 
and is meromorphic in £. By the first resolvent formula, 

R x {K(h n ,0) - Ri(tl(hnXn)) = (l-A)i? A (H(/i„,C))i?i (H(/i„,C)) (28) 
Substituting Eq.(|28|) into Eq.([27). it follows that, 

P(h„, = ^~.i (1 - A)iZ A (H(A n> Q)Ri(tt(h„, 0)dX (29) 
Now, the meromorphic expansion of A i— > Rx(H(h n ,C,)) at Xo(h n ,C,) is, 

i? A (H(/ l „, 0) - t ,, P(h n ,C) + Rx{K, C) (30) 

A - X {h n ,Q 

where A H> R\(h n , () is analytic in the disk {A G € : |A — 1| < 5}. The substitution of Eq.([30]) into Eq.([29]l yields 

P(h n ,C) = (1 - X Q (h n ,())P(h n ,()Ri(m n X)) (31) 
Now, at the pole Cn, the generalized resolvent £ i— > i?i(/i ra ,C) has the meromorphic expansion: 

Ri(h n ,0 = -^—K n + R(h n ,Q (32) 

S Sn 

where H„ is the residue of C n- Ri(h n X) as given in Eq. (fTU|) . and £ i-» R(h n ,() is analytic in the vicinity of Cn- 
Substituting Eg. (|3"2")) into Eq. (f3"Tj) and taking the limit as C — > Cn yields the equation, 

P(/l«, Cn) = -d{X (h n , C„)P(/l n , C«)Hn (33) 



17 



Next, the action of both the sides of this operator equality is evaluated on the state E n . Since P(/i n >Cn) projects on 
E n , and in light of Eq. flTj]) . it follows that, 

E n = — \<p n , E n )d£ \o(h n> Cn)E n 

As the Siegert state E n is not identically zero, it follows that d^Xo(h n , C„) 7^ 0, and the proof of the theorem is 
complete. 

A final noteworthy remark is that the projection Y(h n , £„) and the residue H„ are proportional, which is established 
by applying both the sides of Eq. ([55)l to an arbitrary state F E L 2 (D): 

P(h n ,(n)[F] = -d c \ (h n ,t n )H n [F} 

□ 
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